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Abstract. We consider discrete bilevel optimization problems where the follower solves an integer 
program with a fixed number of variables. Using recent results in parametric integer programming, 
we present polynomial time algorithms for pure and mixed integer bilevel problems. For the mixed 
integer case where the leader's variables are continuous, our algorithm also detects whether the 
infimum cost fails to be attained, a difficulty that has been identified but not directly addressed in 
the literature. In this case it yields a "better than fully polynomial time" approximation scheme 
with running time polynomial in the logarithm of the absolute precision. For the pure integer case 
where the leader's variables are integer, and hence optimal solutions are guaranteed to exist, we 
present an algorithm which runs in polynomial time when the total number of variables is fixed. 

Keywords, bilevel mixed integer linear programming, parametric integer linear programming, 
computational complexity, binary search 

1. Introduction 

Bilevel integer programs are two-stage decision problems in which two decision makers act in 
sequence and where some or all of the decision variables are integer. A leader first decides on 
some variables, denoted here by z (called the leader's, or upper-level, variables), which influence 
the decision of a follower, who then chooses x (the follower's, or lower-level, variables). The leader 
is interested in optimizing over the joint decision (x, z) subject to joint resource and incentive 
compatibility constraints. This setting has been used to model a large variety of applied problems, 
including toll setting in transportation networks [18], revenue management [5], competitive location 
planning [12], and vehicle routing [20] among many others. 

To date the majority of literature has focused on the setting where all constraints and objectives 
are linear and both the leader and follower's variables are continuous (this is evident in observing, for 
instance, the survey papers [3, 4], research monograph [7] and journal special issue [22]). Significant 
work has also been undertaken in nonlinear continuous settings (see for instance [1, 21]). Settings 
where some or all variables are integer are less well-studied, but have nonetheless garnered interest 
in recent years [6, 8, 9, 13, 26]. 
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The main focus of this paper is the following bilevel mixed integer linear program (BMILP), 
where the leader's variables z £ M. d are continuous and the follower's variables x € II 1 are integer: 

inf c • x + e • z (1-1) 

x,z 

s.t. Cx + Dz<p, z>0 (1.2) 
x G argmin x ,{ V ■ x' : Ax 1 <Bz + u, x'eZ"} (1.3) 

where all data are integer. 

An interpretation of the problem is that the leader chooses z and as a result affects the right-hand 
sides of the lower level (follower's) problem 

min{^ • x' : Ax! < Bz + u, x G Z n }, 

x' 

which is a parametric integer program with right-hand side parameterized by z. This interpretation 
underscores the parametric nature of bilevel integer programs that will be a major focus in this 
paper. 

An important comment on the formulation is the implicit assertion of what is known in the bilevel 
literature as the "optimistic assumption" (see for instance [4]). It is possible for the lower level 
problem to have multiple optimal solutions for a given z. The follower is indifferent amongst these 
alternate optima, but the leader may prefer certain optima over others, and yet can exert no more 
direct influence over the follower's choice. How should the follower decide which of the alternate 
optima to choose? To eliminate this difficulty we make the optimistic assumption: whenever the 
follower faces alternate optima for a given z, the follower chooses an optimum x that best suits 
the leader. Indeed, this is reflected in the formulation (1.1)— (1.3) itself: the leader is modeled as 
choosing x within the argmin set of the lower level problem. 

An important object studied throughout is the set of bilevel feasible solutions 

J={(x,z)eK» +i :(1.2), (1.3)}. (1.4) 

We assume that this bilevel feasible set T is bounded (though possibly empty). Much of the 
complication associated with bilevel integer programming is captured in the properties of this set. 
Firstly, T is known to be nonconvex in general, and difficult to describe directly [4]. Secondly, 
given z is continuous it is possible that T is not closed [26]. This possibility is reflected in our 
definition of BMILP in (1.1)— (1.3) where the leader's optimization is expressed as an infimum and 
not a minimum. When z is also integer and the instance is feasible, then an optimal solution always 
exists and "inf" can be replaced with "min" [26]. 

Example 1.1. To demonstrate this complication, consider the following simple example, 
inf —x + z 

x,z 

s.t. 0<z<l 

x G argmin^/j x : x' > z, < x' < 1, i'gZ} 

The "argmin" then simply requires that x = \z\ and thus the problem reduces to: 

infjz - \z\ : < z < 1}. 

2 

Clearly, the infimum is —1 and cannot be attained. Observe that the set of bilevel feasible solutions 

^={(0,0)}lj({l}x(0,l]) 

is neither convex nor closed. 
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The approaches to problem (1.1)— (1.3) in the literature use traditional ideas from integer opti- 
mization, such as branch-and-bound enumeration [23] or cutting planes [9], all based on continuous 
relaxations. It is well documented that there are numerous pitfalls when applying standard ideas of 
branch-and-bound or cutting plane methods to a bilevel setting (see [9] for a discussion). Moreover, 
the algorithm introduced by Moore and Bard in [23] applies only to the case that both x and z 
are integer, or otherwise where it is known that an optimum solution exists. Thus, in general, 
this method cannot solve (1.1)-(1.3), in the sense that it cannot distinguish when an infimum is 
attained or not. The algorithm of DeNegre and Ralphs [9] deals exclusively with the pure integer 
setting. In addition, these two papers give no run-time guarantees. 

Our paper makes the following contributions: 

(i) We present (Theorem 3.1) an algorithm which solves BMILP (1.1)— (1.3) in full generality 
(under our boundedness and integral data assumptions) by identifying when the infimum is 
attained and finding an optimal solution if one exists. No previously proposed algorithm is 
capable of this. 

(ii) Our algorithms runs in polynomial time when the follower's dimension n is fixed for the mixed 
integer case (i.e., z is continuous), and when the total dimension n + d is fixed in the pure 
integer case (i.e., when z is integer). To the authors' knowledge, these are the first polynomial 
run-time guarantees of this type for discrete bilevel optimization. These results can be seen as 
a common extension of two important results on optimization in fixed dimension. The first is 
due to Lenstra [19], who establishes a polynomial time algorithm for integer programming with 
a fixed number of variables. The second is due to Deng [10], who develops a polynomial time 
algorithm for bilevel linear programs when the number of decision variables of the follower is 
fixed. Our result extends Lenstra's result to a bilevel setting and Deng's result to an integer 
setting. For further discussion of complexity issues in bilevel optimization see [3] . 

(iii) When no optimal solution to BMILP exists, we find, for any given e > and in time polynomial 
in the input size and log(l/e), a solution with objective value no more than 1 + e times the 
infimum. 

In Section 2, we introduce formal definitions and prove our first key result: when an instance 
of BMILP is feasible its infimum value is a rational number with polynomially bounded binary 
encoding size. In Section 3, we apply recent results from the study of parametric integer program- 
ming by Eisenbrand and Shmonin [11], based on the work of Kannan [14]. These references provide 
polynomial time algorithms for decision versions of parametric integer programs. We apply these 
algorithms in the mixed integer bilevel setting and, using binary search, derive our main result 
(Theorem 3.1): a polynomial time algorithm for the mixed integer case where z is continuous and 
the follower's dimension n is fixed. Finally, in Section 4 we present two algorithms to find optimal 
solutions to the pure integer case; that is, BMILP under the further restriction that z is integer. 
Both of these algorithms run in polynomial time when the total dimension n + d is fixed. The first 
algorithm is similar to those discussed in the previous section, whereas the second is based on an 
approach similar to that found in [15] which is based on a theory of rational generating functions. 



2. Definition and Preliminary Results 

An instance of the bilevel mixed integer linear program (BMILP) (1.1)-(1.3) is defined by the 
leader's and follower's dimensions d and n; the numbers h and m of linear inequality constraints 
in (1.2) and in the feasible set of the follower's subproblem (1.3); integer matrices A G ^rnxn^ 
B G Z mxd , C G Z hxn and D G Z hxd ; integer objective coefficient vectors c G Z n , e G Z d and 

G Z"; and integer right-hand side vectors u G Z m and p £ Z*. The binary encoding input size 
of an instance is the total number of bits needed to input all the data defining the instance. 

3 



We assume that the upper level polyhedron 

P = {(x,z) GR n+d : (1.2)} (2.1) 

is bounded (but possibly empty). We also assume that for every z G proj z P, where 

proj z P = {z G R d : (x, z) € P for some x G M n }, 

the lower level feasible set {(x, z) G W n+d : Ax < Bz + u} is also bounded (and possibly empty). 
(For example, by Farkas's Lemma, if lower bounds x > x and z > z are part of or implied by the 
constraints defining P, and also by those defining the lower level feasible set, then the boundedness 
assumptions are satisfied if and only if the (n + <i)-row vector of all l's is a nonnegative combination 
of the rows of the constraint matrix (C, D), and also a nonnegative combination of the rows of the 
matrix (A, —B).) 

Under these boundedness assumptions, an instance of BMILP is either infeasible or else must 
have a bounded infimum value. Feasibility can be decided in polynomial time by applying Lenstra's 
(mixed) integer programming algorithm [19] to decide if {(x, z) G P : Ax < Bz + u, x G Z™} ^ 0. 
A key result used in the sequel is that, when the instance of BMILP is feasible, its infimum value 
is a rational number with polynomially-bounded binary encoding size. 

Proposition 2.1. Assume that a given instance of problem (1.1)— (1.3) is feasible and has a bounded 
optimum (infimum) value. Then the infimum is a rational number with denominator bounded above 
by the maximum absolute value of a sub -determinant of the matrix (D,B) T . 

Proof. Let (x, z) G T be a bilevel feasible solution to (1.1)— (1.3) and define 

Q(x, z) = {zeR d :Dz<p- Cx; #(z) < B { z + m < #(z) + 1 Vi G Af }, 

where M = {1, . . . , m}, B{ is the i-th row of B, and /3j(z) = \_BiZ + m\ . Note that z G Q(x, z), so 
Q(x, z) is a nonempty, bounded quasi-polyhedral set, i.e., its topological closure is 

clQ(x,z) = {z G R d : Dz < p - Cx; &(z) < Bp + Ui < ft(z) + 1 Vz G M}, (2.2) 

a nonempty polytope. Every z G Q(x, z) satisfies \BiZ + UjJ = \ Biz + tijj for all i = 1, . . . , m and, 
since A is integer, 

VxGZ": Ax < Bz + u Ax < Bz + u, 
implying that x G argmin{-0 • x : Ax < Bz + u, x G Z n } for all z G Q(x, z), that is, 

VzGQ(x,z) (x,z)GJ zr . (2.3) 

Now let v* denote the infimum value in BMILP (1.1)— (1.3). There is a sequence (x fc , z k ) of 
bilevel feasible solutions with objective value c • x + e • z converging to v*. Since P is a polytope, 
its projection proj x P = {x G 1" : (x, z) G P for some z G M. d } contains a finite number of integer 
points x G -PnZ n . Therefore there is a subsequence (x fc ,z fc )fc e N with all x fc = x°. For every i 6 M 
there is also a finite number of possible integer values j3%{z) = \_B%z + Ui\ when z G proj z P. Thus 
there exists a subsequence (x , z )fcgN m 

^ with all (3{z k ) = P(z°), i.e., (x°,z fc ) G Q(x°,z°), and 
c • x° + e • z k converging to v* . From (2.3) it follows that v* = c ■ x° + inf z {e • z : z G Q(x°, z )}. 
By continuity of the objective e • z, v* — c ■ x° is the optimum value of the linear programming 
problem min z {e • z : z G cl Q(x°, z )}. By (2.2), for every i G M at most one of the two constraints 
— -^« z + u i or + n « — A(z°) + 1 is binding at the LP optimum, hence v* — c • x° must 
also be the optimum value of a linear programming problem 

min e • z 

s. t. Dz < p - Cx° 

- B iZ < Ui - (3i(z°), ViG/ (2 ' 4) 
Pjz < - Uj + /3j(z°) + 1, Vj G J 
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for some disjoint I, J Q M. It is well known that the optimal value of a bounded linear program 
with integer data occurs at a basic feasible solution, the denominators of which are bounded by 
the maximum absolute value of a sub-determinant of its constraint matrix. By simple properties 
of determinants the stated result then follows. □ 

Note that in the pure integer case, where the leader's variables are also required to be integer, 
the bilevel feasible set is 

T 1 = {(x,z) G P : z G Z d , x G argmin{ 1/7 • x' : Ax! < Bz + u, x G Z n }} . (2.5) 

When it is non-empty, T' is a finite set, and then the infimum in (1.1) is attained. Since c and e 
are integral, this optimum value is integral, providing a trivial strengthening of Proposition 2.1 in 
the pure integer case. 

3. An Algorithm for Bilevel Mixed Integer Linear Programming 

In this section we prove the main result of this paper: 

Theorem 3.1. There exists an algorithm which solves the BMILP problem (1.1) -(1.3) in the 
following sense: the algorithm 

(i) decides if the instance is feasible; 

(ii) if it is feasible, decides if an optimal solution exists (that is, if the infimum is attained); and 
(Hi) if the infimum is attained, finds an optimal solution. 

The algorithm runs in polynomial time when the follower's dimension n is fixed. 

To prove this theorem, we use results on parametric integer programs with parameterized right- 
hand sides. Parametric integer programming is well-studied, the particular approach we discuss 
builds on geometric and algorithmic number theory and has it roots in the work of Lenstra [19] on 
integer linear programming in fixed dimension. We draw mainly on the recent work of Eisenbrand 
and Shmonin [11], which itself builds on earlier results by Kannan [14]. These two related works 
solve a decision version of parametric integer programming (defined below) in polynomial time 
when certain parameters are fixed. As discussed in Section 1, we can view the lower level problem 
as a parametric integer program. But the connection between parametric integer programming and 
bilevel integer programming in fact goes deeper, as we now discuss. 

The integer projection of a polyhedron Q C M m+P is the set 

Q/Z p = {b£R m : (b, w) G Q for some w G IF}. (3.1) 
The feasibility version of parametric integer linear programming is the following decision problem: 
PILP: Given a rational matrix R G Q mxn and a rational polyhedron Q G M m+P 
decide if the linear system Ry < b has an integer solution for all b G Q/1F. 

Letting 

Pb = {y € R n : Ry < b}, 
this problem is equivalent, by negation, to deciding the sentence 

3heQ/Z p p b nz n = 0. (3.2) 

Eisenbrand and Shmonin [11] prove: 

Theorem 3.2. (Theorem 4.2 in [11]) There is an algorithm that, given a rational matrix R G Q mxn 
and a rational polyhedron Q C ]R m+p , decides the sentence (3.2). The algorithm runs in polynomial 
time if p and n are fixed. 

Our goal is to establish a connection between PILP and BMILP so this result can be used in 
solving BMILP. First we define a decision version as follows: 
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BMILP a : Given integer matrices A G Z mxn , B G Z mxd , C G Z hxn and Z) G Z hxd ; 
integer vectors c G Z n , e G Z d , ip G Z n , u G Z m and p G Z ft ; and rational scalar 
a G Q, decide if there exists a vector (x, z) G M n+d such that (1.2), (1.3) and 
c • x + e • z < a. 

Proposition 3.3. There exists an algorithm that decides BMILP Q . The algorithm runs in polyno- 
mial time when the follower's dimension n is fixed. 

Proof. Since ijy is integral, we can express (1.3) as the fact that 

{x'£r:f x'<f x-1, Ax' < Bz + u} n Z n = 0. (3.3) 

Let R = (-0, A, Odxn) T where Odxn is a d x n matrix of zeroes; b = (&m, t>(2)> b(3)) T with bm = 
■0 • x — 1, b(2) = Bz + u and bp) = z. Then (1.3) and z > together are equivalent to P^ n Z n = 
where P£ = {x' G W 1 : i?x' < b}. Letting 

Q Q = {(b,x) : c • x + e • b (3) < a, Cx + Db (3) < p, b (3) > 0, 

v4x<b (2 ), 6(i) =-0 -x-1, b (2 ) = Bb( 3 ) + u}, (3.4) 

deciding BMILP a reduces to deciding an instance of (3.2) with p = n, w = x, Q = Q a and 
P\j = Pi. The result now follows from Theorem 3.2. □ 

Our approach to the proof of Theorem 3.1 is to use binary searches on the target value a in 
BMILP a and related decision problems, so as to find the infimum value of problem BMILP and an 
optimal solution if one exists. 

The following result from Kwek and Mehlhorn [17] yields the binary search algorithm used in 
proving our results: 

Theorem 3.4. (Theorem 1 in [17]) There exists a 9(\ogL)-time algorithm that determines an 
unknown rational number r which lies in the set Ql = || : p, q G {1, . . . , £}| by asking at most 
2 log 2 L + 0(1) queries of the form "is r < a ?". 
We first establish: 

Lemma 3.5. Consider a feasible instance to BMILP problem (1.1) —(1.3) . Its infimum value can 
be determined in polynomial time when the dimension n is fixed. 

Proof. We first determine polynomial-sized lower and upper bounds v and v on the infimum value 
v* by solving in polynomial time the linear programs v = min{c • x + e • z : (x, z) G P} and 
v = max{c • x + e • z : (x, z) G P}. We then perform ordinary binary search for v* in successive 
intervals [u',u"] satisfying v' < v* < v", where the search query "is v* < a" is the decision problem 
BMILP Q , and the initial interval [v' , v"\ = [v — l,v]. We stop, after a polynomial number of queries, 
as soon as v" — v' < 1. If there exists an integer j3 such that v' < (3 < v" then we perform one 
more query "is v* < (5" and set the current interval to \J3 — l,/3] if the response was positive (i.e., 
if v* < (3), and to v"\ otherwise. If there is no such integer (3 then we set the current interval 
to [Lf'J,f"]. In any case we obtain an interval [u', - ;;"] such that v' < v* < v" , v" — v' < 1 and 
v' is integer. We now invoke the binary search algorithm of Theorem 3.4 for the rational number 
r = v* — v' G Ql where L is the upper bound on the denominator of v* from Proposition 2.1. □ 

We are now ready to prove Theorem 3.1. 

Proof of Theorem 3.1. (i) As indicated in Section 2, feasibility can be decided in polynomial time 
by applying Lenstra's (mixed) integer programming algorithm [19] to decide if {(x, z) G P : Ax < 
Bz + u, x G Z n } ^ 0. 

(ii) Assume the instance is feasible. By Lemma 3.5 we can find the infimum value v* in polynomial 
time. This infimum is attained if and only if there exists a bilevel feasible (x, z) G T such that 
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c • x + e • z = v* . This can be decided in polynomial time as in the proof of Proposition 3.3 by 
replacing the first constraint in the definition of Q a in equation (3.4) with c • x + e • b( 3 ) = v*. 

(iii) Finally, assume that the infimum value is attained. We extend to bilevel mixed integer 
programming a standard method used in (ordinary) mixed integer programming to successively 
determine an optimum solution given an optimum objective value oracle (see for instance [24], page 
260). In this method, we first find the optimum value of related instances of BMILP to successively 
determine each component of an optimal solution (x*,z*) with lexicographically smallest integer 
component x*. By induction, assume that all components x* k for k = 1, . . . , j — 1 have already been 
determined. Then x* is the optimum value of the BMILP 

min Xj 
x,z J 

s.t. Cx + Dz < p; z > 0; c-x + e-z = t>*; 
x k = x% for k = 1, . . . ,j - 1; (1.3). 

Since this problem is feasible and its objective is integer valued and bounded, its optimum value is 
attained, and is determined in polynomial time by Lemma 3.5. (Note that, since this optimal value 
x* is integer, we only need to use ordinary binary search and stop when the width of the current 
interval is less than one.) 

Having determined x*, we next determine the continuous component z*. Note that, even though 
x* has now been fixed, z* cannot be directly found by solving a simple linear program, because 
we still need to enforce the (nonconvex) incentive compatibility constraint (1.3). Using ideas from 
the proof of Proposition 2.1, we construct an "integer RHS vector" r = (LBjZ + UiJ)i=i,... m which 
is lexicographically minimum among all those defined by optimum BMILP solutions (x*,z). Thus 
by induction assume that for k = 1, . . . , i — 1 we have determined integers such that the set 

R _ J '/v v \ <= i$n+d , Cx* + Dz<p; z>0; c-x* + e-z = t;*; 1 
<j^x,z;tK . ^ < BkZ + Uk <rk + i for all fc < i; x = x*; (1.3)/ 

is nonempty. As in the proof of Proposition 2.1, note that the closure cl-Rj_i is determined by 
replacing each strict inequality B^z + < + 1 with B^z + < ru + 1, and is thus a nonempty 
polytope. We then compute the infimum value of the BMILP 

Pi = mf {BiZ + Ui : z £ dRi-i}, (3.5) 

x,z 

in polynomial time by Lemma 3.5. Let T{ = [pi\ . We claim that Ri ^ 0. To prove this claim, let 
((x*,z J ) G R%-i) gN be a sequence of points in cli2j_i with limj_ +DO (.Bj z? +Uj) = p^. By continuity 

of the linear function (x*,z) i— > Bi z + itj, there exists a sequence ((x*,z J ) G Ri-i) - eN of points 
in Ri-i with lim J _ +00 (Bj z J + Uj) = pi < [pi\ + 1. This implies that |_PiJ < A z J + u- L < [pi\ + 1 for 
some j G N, and therefore (x*, z-?) G Ri, showing that Ri ^ 0. 

Thus, after determining x* and computing the infimum value pi of m BMILPs (3.5), we obtain 
the polynomial-sized, integer vector r = m such that there exists an optimal solution 

(x*,z) to the original BMILP (1.1)-(1.3) satisfying 

. „ _ „ f md Cx* + Dz < p, z > 0, c • x* + e • z = v* 1 
zGpro Jz i? m CQ:=|zGM d : r < ^ + + f ) 

where 1 is the all-ones vector in Z m . As in the proof of Proposition 2.1 we note that every z G Q 
gives rise to an optimal solution (x*,z). We now show how to find one such optimal solution 
whose binary encoding length is polynomial in the input size. Let k = 1 + dimQ. Find k affinely 
independent vertices z 1 , . . . , z k of cl Q. This can be done in polynomial time since cl Q is a polytope 
(see, e.g., [24], page 186). By standard results in linear optimization, the vertices z 1 , . . . ,z k can be 
written as rational numbers in ^L d where A is the least common multiple of the determinants of 
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the bases of the constraint submatrix associated with the basic feasible solutions z , . . . , z . Note 
that A has encoding length polynomial in the input size. It is possible that none of the i? lie in 
Q, since Q is not closed, but nonetheless their barycenter z* = Y2j=i \ ^ ls guaranteed to be in Q 
and also in the set -^^ d - Thus it has a binary encoding length polynomial in the input size. 

It follows that (x*,z*) is an optimal solution to (1.1)-(1.3) and can be found in polynomial 
time. □ 

We note, in particular, that this theorem provides a polynomial-time algorithm to decide if a 
BMILP instance has an optimal solution. To the authors' knowledge, this is the first algorithm to 
achieve this in the literature on discrete bilevel optimization. The existence of solutions in the case 
where z is continuous has long been known not to be guaranteed (see for instance [23]) but little is 
known in terms of necessary and sufficient conditions for existence. Dempe and Richter [8] explore 
the setting where the lower level problem is a knapsack problem, and they provide necessary and 
sufficient conditions for existence of solutions, but there have been no studies of this kind in a 
more general setting. An efficient algorithm to decide the existence of solutions, such as the one in 
Theorem 3.1, may be helpful in investigations of this kind. 

A final question remains: what can we say about the BMILP problem (1.1)— (1.3) when the 
infimum is not attained? In such a case it may be desirable to find e-optimal solutions whose 
objective values approximate the infimum. The following result demonstrates how we can adjust 
our arguments to find e-optimal solutions. 

Corollary 3.6. Let e > be given in binary encoding. Given a feasible instance of the BMILP 
problem (1.1)— (1.3) with infimum value v* , there exists an algorithm to find a feasible solution (x,z) 
that satisfies c-x + e- z<(l + e)v* . The algorithm runs in polynomial time when the follower's 
dimension n is fixed. 

Proof. By Theorem 3.1 part (ii) we can decide in polynomial time if the infimum is attained. If 
so, by using part (iii) of the same theorem we can find an optimal solution (x*,z*) that satisfies 
c-x* + e-z* = v* < (l + e)v*. 

If we find that the infimum is not attained we apply a similar procedure as in the proof of 
Theorem 3.1 part (iii). By replacing the condition c • x* + e • z = v* with the condition c • x + e • z < 
(1 + e)v*, the result follows immediately. □ 

It should be noted that this algorithm achieves a stronger form of efficiency than a fully polyno- 
mial time approximation scheme (FPTAS). In an FPTAS, the error e is given in unary encoding, 
whereas in our algorithm e enters as input in binary encoding, since its value only plays a role in the 
constraint c-x + e- z < (1 + e)v*+. In this sense, we have obtained a "better than fully polynomial 
time" approximation scheme to yield solutions with objective value e-close to the infimum. 

A distinguishing feature of our main result, Theorem 3.1, is that the dimension d of the leader's 
variable need not be fixed in order to assure a polynomial running time. Note, however, that this 
result in no way conflicts with the fact that bilevel linear programming, in which both x and z 
are continuous and the dimension of both x and z are allowed to vary, is an NP-hard problem. 
In our setting it is the leader's variable only which is allowed to be continuous and have varying 
dimension, and an important ingredient in our proof is that there are only finitely many integer 
right-hand sides /3j(z) = [B^z + uj to consider when the matrix A and decision variables x are 
integer. 

4. An Algorithm for Bilevel Integer Linear Programming 

We can also extend our results to the pure integer setting, where z is further restricted to be 
integer. We refer to this problem as the bilevel integer linear programming (BILP) problem, an 

8 



instance of which is specified by the same data as BMILP and the same constraints (1.1)— (1.3) with 
the additional integrality condition 

z G Z d . (4.1) 
As above, we assume that the feasible set T' defined in (2.5) is bounded. 

Theorem 4.1. There exists an algorithm that solves the BILP problem (1.1)-(1.3) & (4.1) in the 
following sense: the algorithm 

(i) decides if a given instance is feasible; 

(ii) if it is feasible, finds an optimal solution. 

The algorithm runs in polynomial time when the total dimension d + n is fixed. 

Proof, (i) Feasibility can be established in polynomial time by applying Lenstra's integer program- 
ming algorithm [19] to decide if {(x, z) € P : Ax < Bz + u, (x, z) G Z n+d } ^ 0. 

(ii) Assume the instance is feasible. The idea of how to determine an optimal solution is similar 
to, but simpler than, the mixed integer case, as we use binary search to directly find the lexi- 
cographically minimal optimal solution (x*,z*). To begin, consider the decision problem BILP Q , 
which is identical to BMILP a with the additional integrality condition (4.1). 

Problem BILP a can be decided in polynomial time using an algorithm similar to that in Propo- 
sition 3.3. Indeed, if we simply redefine 

Q a = {(b,x,b) : c • x + e • b (3 ) < a, Cx + Db (3 ) < p, b (3) > 0, 
Ax<b (2 ), b (1 )=y-x-l, b (2 ) = 5b( 3) + u}, 

and let p = n + d, the algorithm in the proof of Proposition 3.3 decides BILP a . 

It is then straightforward to find the optimal value v* of the given BILP instance using binary 
search as in the proof of Lemma 3.5 and stopping when the search interval has length less than 
one. 

Given the optimal value v* we again use binary search to determine the lexicographically minimal 
optimal solution (x*,z*). The procedure to find x* is identical to that described in the proof of 
Theorem 3.1. The procedure to find z* is also similar. We successively determine each component 
of z* as follows: assume that all components z k for k = 1, ... ,j — 1 have already been determined. 
Then Zj is the optimal value of the BILP 

min Zj 

z J 

s.t. Cx* + Dz<p; z > 0; z € Z d c ■ x* + e • z = v*; 

z k = z* k for k = 1, - 1; 

x* G argmin x ,{ V ■ x' : Ax.' < Bz + u, x' G 7L n }. 

Since this problem is feasible and its objective is integer valued and bounded, its optimum value is 
attained, and can be determined as in the previous paragraph. 

The result then follows. □ 

We point out that an alternate algorithm satisfying the conditions of the previous theorem can 
be obtained by applying a rational generating function technique similar to that found in [15] . The 
algorithm is very similar to that in the Stackelberg-Nash setting of Section 6 in that paper. The 
differences are: (i) there is only a single follower as opposed to a group of followers (hence any 
follower's optimal decision is trivially a "Nash equilibrium"); (ii) the payoffs are linear instead of 
piecewise linear; and (hi) the set constraining the leader's actions is more complicated and depends 
on the action of the follower, see (1.2), whereas in the Stackelberg-Nash setting of [15] the leader's 
choice set is independent of the actions of the followers. Nonetheless, these differences are minor 
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and easily handled by the rational generating function technique described in [15]. Details are thus 
omitted. 

How do these two algorithms to solve BILP - the one based on parametric integer programming 
versus the one based on rational generating function techniques - compare in terms of efficiency? 
It is accurate to say that the algorithm using parametric integer programming programming is 
more direct and simpler. Indeed, the rational generating function approach relies on a strong result 
called the "Projection Theorem" (Theorem 1.7 in [2]) which, although it has polynomial running 
time (in fixed dimension), was until recently thought to be virtually unimplementable in practice. 
Only recently has the Projection Theorem been implemented in the library barvinok developed by 
Sven Verdoolaege [25]. Moreover, the current implementation appeals to the results in [11], some 
of the very same results that we use more directly in the development of our first approach (see [16] 
for more details). Future work would be required to improve implementations of the Projection 
Theorem to yield a practical solution method for bilevel integer programs. 

As for computational comparisons, in the pure integer setting, between the our two algorithms 
of Theorem 4.1 herein and those of Moore and Bard [23] and DeNegre and Ralphs [9], this is also 
a direction for future work. 

5. Conclusions 

In this paper we have detailed an algorithm which solves the BMILP problem (1.1)— (1.3) in 
polynomial time when the dimension of the follower's integer decision variable is fixed. It solves 
BMILP in the following strong sense: under our boundedness and integral data assumptions, it 
decides if the problem is feasible, and if so determines its infimum; if the infimum is attained, it 
finds an optimal solution, else it finds an e-optimal solution for any e > given as input in binary 
encoding. 

This algorithm is based on recent results in parametric integer programming due to Eisenbrand 
and Shmonin [11] combined with binary search to find the infimum objective value and an optimal 
or e-optimal solution. A simplified version of our algorithm also solves the pure integer problem 
BILP (1.1)-(1.3) & (4.1) in polynomial time when the total number of decision variables is fixed. 
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